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ABSTRACT 


Spaceborne altimetry with measurement accuracies of a few centimeters has the 
potential for determining s:'a-surface elevations necessary to compute accurate three- 
dimensional geostrophlc currents from traditional hydrjgraphlc observations. The 
limitation In this approach Is the uncertainties In cur knowledge of the global and 
ocean geopotentials which produce satellite and gecid height uncertainties about an 
order of magnitude larger than the goal of about ’0 cm. This paper begins with a 
description of the quantatlve effects of geopotentlal uncertainties on processing 
altimetry data. This is followed by a review of existing modtia wmd. shown to 
be Inadequate. Potential near-term Improvements, not requiring additional spacecrau, 
are discussed. However, even though there would be substantial Improvements at the 
longer wavelengths, the oceanographic goa 1 would not be achieved. The potential NASA 
Geopotential Research Mission (GRM) Is described. This mission should produce geo- 
potential models that are capable of defining the ocean geoid to 10 cm and near- 
earth satellite positions significantly better. For completeness, the state-of-the- 
art and the potential of spacebornr gravity gravimetry is described as an alterna- • 

tlve approach to Improve our knowledge of the geopotential. j 

! 

1 . INTRODUCTION 

* 

This paper addresses the Importance of an accurate representation < * the geo- 1 

potential In physical oceanography. Its current state oi knowledge and possible near- i ! 

term and long-term Improvements. : 

; t 

If the oceans *»crs static and subject only to gravitational and centrifugal 1 

forces, they would conform to a conceptual surface called the geoid on which the ! ; 

gravity potential function (or geopotentlal ) is a constant. The oceans are not 
static, and transport of heat, salt, and morentum have a profound effect on climate. 

As a result, study of the temperal and ste-jdy-state circulation of the oceans Is • 

Important. These motions are governed bv the equations of fluid dynamics and solu- A - 

tlon Is dependent on conditions at the boundaries and throughout the medium. 

Velocities at the ocean surface are primarily Induced by wind stress *nd differences ! 

In pressure. On the rotating earth, ocean pressure differences give wise to geo- • * 

strophlc velocities resulting from a balance between pressure and corlolls forces ( 

such that the velocity Is normal to the pressure gradient. Surface pressure differ- 
ences are caused by departures of the ocean surface from the ocean geoid. These 
departures, called sea-surface elevations, can be as large as 1 to 2 meters In the j 

broad ocean areas and significantly larger at the land-sea Interfaces. \ 

The relation between the surface geostrophlc velocity and change In s ^-surface • 

elevation Is given by j 

v s (g/f)(3H/3x) ‘ j 

where g Is the local acceleration of gravity (*9.8 m.s" 2 '-, ? ■ 2o sin (latitude) where • , 

a Is the angular velocity of the earth (7.27 x 10" 5 £'*), and 1' che horizontal 
slope of the sea-surface elevation. An uncertainty In sea- Si*rf -.*.»• elevation of 10 cm 
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over a distance of i 000 km would result in a velocity uncertainty of 1 cm-s* 1 at 
a latitude of 45 degrees. This would be sufficient for a better understanding of 
• current velocities. 

A gross quantitative understanding of ’arge-scale ocean circulation has resulted 
from shipboard measurements. Velocities based on density distributions inferred 
from in situ measurements can be determined to within a local constant if the sea- 
surface elevations are unknown. Determination of this constant is possible by in situ 
measurement of velocity. However, in practice this is not usually feasible as it 
would take months to average the small-scale velocity variations. Alternatively, the 
• constant can be determined from the suspect assumption that the velocity vanishes 
at some level in the ocean called the level-of-no-motion. The inability to accurately 
determine this constant has been a limitation to the quantitative description o t ' the 
geostrophlc current systems. 

With the advent of spaceborne microwave altimet r y a solution is possible. With 
its inherent global coverage and measurement accuracy, altimetry has the potential 
for determining sea-surface elevations to u few centimeters. This, together with tra- 
ditional hydrographic observations, would make possible determination of the three- 
dimensional geostropnic currents unencumbered by the assumption of a level-of-no-motion. 
However, the effectiveness of this approach Is limited by uncertainties in our know- 
ledge of the global geopotential and the ocean geoid. 


2. SATELLITE ALTIMETRY 

Spaceborne nadir pointing high resolution microwave radar altimetry has been an 
exciting source of data for ocean topography (Chovitz, 1983; and Marsh, 1983). 
Instrument range measurement precision has improved steadily from the 1 to 2 m for 
Skylab and the 30 to 40 cm for GECS-3 to the 5 to 7 cm for SEASAT. Two new space- 
borne altimeter missions are planned. The U.S. Navy program, GEOSAT, is undergoing 
fabrication at The Johns Hopkins University Applied Physics Laboratory with launch 
scheduled in the fall of 1984 (Pisacane and DeBra, 1983). Primary purpose of the 
mission is to better determine the ocean geoid by essentially completing the SEASAT 
mission. By collecting data at widely spaced Intervals of time errors caused by 
time dependent sea-surface elevations, can be minimized. Instrument accuracy 
should be the same as SEASAT because it is essentially of the same design with im- 
provements of an engineering nature. The second program Is TOPEX which is a NASA 
venture now in the planning and instrument development stage. The primary purpose 
of TOPEX is to determine sea-surface elevations. An official new start is projected 
for 1985 with launch expected in 1988 or 1989. 

Figure 1 shows the geometrical configuration of satellite altimetry. From 
the altimetry data, h can be Inferred where K is the distance from the center of 
mass of th$ spacecraft to that point on the ocean nearest to it. Sea- surface 
elevation H is given in terms of n by 

H*r-fi-N-r g 

where N is geoid height vector and r and r g are the position vectors to the center of 

* . mass of the spacecraft and to the subsatellite point, on the reference surface respec- 
* tlvely. Errors in «H are 

6H * {j-(6r - fit - 6N) 


*1 
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Position vectors from centjr of mass reference frame: 
T = satellite 

Fg * reference spheroid 

Height vectors: 

h = altimeter measurement 
fl * sea surface elevation 
N * geoidal height 

Fig. 1 Satellite altimetry geometry. 


Because the vectors are near parallel. 


fiH = fir - fih - fiN 

so that errors In the satellite altitude, the measurements, and the geold height have 
the same sensitivity. Because of the large footprint of the altimeter, i.e. , v 1 km, 
errors In spacecraft position orthogonal to n are of second order. Consequently, r, 
h, and N all need be determined to the same degree of accuracy. 
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Contributions to 6h are uncertainties in instrumentation delays, distance from 
the center of mass to the electrical center of the antenna, spacecraft attitude, 
propagation velocity, the effect of ocean surface characteristics, and random noise. 
These are discussed in depth by Marsh (1983) and Tapley et al (1982) and will not be 
considered here as ether errors dominate. Errors that contribute to fir, i.e., satel- 
lite altitude, are of two types: those from the spacecraft tracking system (e.g., 

station location, propagation velocity, and instrumentation errors) and modeling 
of the forces that act on the spacecraft, (e.g., gravity, radiation pressure, and 
drag) that are necessary to correlate tracking data taken at different times. Today, 
the dominant error in determining satellite ephemerides from tracking systems such 
as laser and radiofrequency doppler is the uncertainty in the global geopotential. 
This uncertainty manifests itself in tracking station position errors and errors in 
the gravity forces. Errors in the geoidal height, 5N, follow directly from both 
errors in and truncation of models of the geopotential. Because of density inhomo- 
geneities, primarily in the earth's crust, the geoid is not a smooth surface and can 
depart from the reference oblate spheroid surface by as much as 100 m. 


The total gravitational potential, V*, can be represented in terms of spherical 
harmonics by 

V*(r,A,*) = V 0 + V, 


where 


r m 

= — , the Newtonian potential; 


( 1 ) 


£m 


= E 


$ v 


1=2 m =0 

5M ,R 
R l r 


Am 


GM (?) £+1 (C._ Cos nu + S 0m Sin mx) P„ m (cos *) 
v £m £m £rrr 


( 2 ) 


G 

M 

R 


= universal gravitational constant 


= mass of earth 


* normalizing radius, generally the mean equatorial radius 


= radius, longitude and colatitude 


Am 


= (1 - t^) m ^ -~P (t), associated Legendre function of the 


dt' n 1 first kind 


P # (t) a — ^7 (t^ - ]) 1 , Legendre polynomials 

1 2 l i\ dt* 


'to 


jlj- /(£ ) A P^(cos ♦ * )dM 
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r^A 1 ,^' = integration parameters over the mass of the earth 

CjmT ^£m are t ^ ie harmonic coefficients which are unknown integrals of the mass 

distribution of the earth. Following Dunnell et al (1977), the spherical harmonic 
V can be represented in Kepler elements (a,e,i for near zero eccentricity 

by 


V 


£m 


GM 

~T 
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t+i 


£ 

I I nmri S_ 
p=0 *"ip £mp 
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Where 



r 

(£-m)even 

S £m 

r 

(£-m)even 


S £mp 

L £m 

_c 

cos ♦ + 

£mp 

sin 4> . 

(4) 


£m 

(£-m)odd 

L £m 

(£-m)odd 



%p = + m • e )- ( 5 ) 


a is the semimajor axes, I is a function of the inclination, b ■ M+cu is the argu- 
ment of latitude, n is the longitude of the ascending node, and e Is the right 
ascension of Greenwich relative to Arles. This provides a convenient representation 
of the geopotential for deriving both the geold height and the effect of the harmonic 
coefficients on satellite motion. Geold undulations can be determined from the Bruns 
theorem as 


N 



( 6 ) 


o 

where g s GM/R Is the local acceleration of gravity and U is the difference of the 
potential of the ellipsoidal reference surface with the Newtonian potential for which 
all terms for £ > 2 can be neglected. Substituting for V and approximating a by R 
gives 
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Perturbations In the radial direction of a satellite In a near-circular orbit, 
following Dunnell et al (1977), are 
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( 8 ) 


Because of the linear nature of the expressions for N £ and <$r^ in terms of the 
harmonic coefficients and S £m> equations (7) and (8) can also be interpreted 
as determining the effect of uncertainties in the harmonic coefficients. 

Estimates of the relative values of 6r £ and N £ can be determined as follows. 
If the earth were essentially nnnrotating such that n-e is a constant then 


*tmp ■ 

2 3 

and using n = GM/a , equation (8) can be written as 


«r » R l (7) 
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The maximum contribution to the perturbation from the harmonic coefficients at a 
given frequency + (£-2p)n occurs when £ is a minimum for which p must be either 0 
or £. For the highest frequency in -Sr, denoted by 6r ( + , the perturbation is 


6r £|*=+£n = R( a )4 " 1 J n FT [I £mo s wno + WW 

1 — FTFO 


(10) 


Similarly, for N„ 
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£ *»+in 
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( 11 ) 


Then the amplitude or the uncertainty of the ratio of the perturbation In satellite 
height to the geoid undulation at a frequency £n is 


* r i ; 


*=+£n 


<f> 


£-1 


(A + 1) 


-1 


(12) 
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This ratio as a function of harmonic degree for two satellite altitudes Is 
given In Figure 2. Attenuation of the ratio as £ Increases demonstrates the diffi- 
culty in determining small-scale variations In the geopotential by measurement of 
orbit perturbations. This also demonstrates that the orbit is less effected by 
contributions of the harmonic coefficients at higher degree £. Rotation of the 

such that the 


:arth Introduces frequencies that are smaller than orbital, < n. 


magnitude of the orbital perturbations are enhanced by the factor « 


1- 1 


£mp* 


Consequently, 


the more significant perturbations occur at the longer wavelengths and small-scale 
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altimeter variations are primarily a result of the topography of the ocean surface. 
Equations (10) and (11) demonstrate that small variations in the mean Kepler elements 
produce second-order effects in both the satellite altitude perturbation and the 
ocean geoid. This is the mathematical justification of using an orbit with a re- 
peating ground track so that the geopotential uncertainties result in highly cor- 
related errors. Variations in the altimeter measurements can then be interpreted 
as time dependent components of sea-surface elevations. 


To determine the absolute sea surface elevation requires a global geopotentiai 
model to determine the satellite altitude which is far less detailed than the model 
required to define the ocean geoid. 


Wevelenath 
lower oound 
(km) 
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Fig. 2 Ratio of satellite height error to geoid height error 
versus harmonic degree. 


3. GEOPOTENTIAL MODELING STATUS 

A review of the current status of modeling the geopotential is available (Lerch 
1983). These models use various sources of data: orbit perturbations through 
laser and radio frequency doppler observations} satellite altimetry, and measurements 
of terrestrial gravity anomalies. Resolutions are as small as a half-wa/elength of 
one degree. 

Over the years, NASA has generated the Goddard Earth Model (GEM) series which 
has been the accepted standard. The latest In the series are GEM-9 and 10 (Lerch 
et al, 1979), GEM-10B and IOC (Lerch et al, 1981), and GEM-L2 (Lerch et ai, 1983). 
Characteristics of these models are given in Table 1. GEM-9 and GEM-L2 are based 
solely on satellite tracking data, GEM- 10 Is based on the data used In GEM-9 but 
augmented by terrestrial gravity anomaly measurements, and GEM-10B and 10C use 
satellite tracking, surface gravimetry and GEOS-3 altimetry. Various measures of 
accuracy for GEM-9, 10, 10B and 10C are given In Table 2. 
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For orbit determination GEM-9, 10, 1 OB and 10C are comparable with a radial 
accuracy of about 1 to 2 meters. This compares favorably to the 1 to 2 meters esti- 
mated for the satellites of the Navy Navigation Sa'ellite System at about 1000 km. 

The GEM-L2 model is reported to be superior to the i arlier models by more accurate 
determination of the harmonic coefficients up to degree and order 4. The orbit 
errors for the LAGEOS spacecraft at an altitude of almost 5900 km were reduced 
from 2m with GEM-9 to 0.30 cm for GEM-L2. There has not yet been evaluation of 
this model for low altitude satellites. 

The measures of accuracy given in Table 2 indicate that the GEM-10C model has 
an edge in reproducing the ocean geoid. An accuracy of 1 to 2 meters in geoid height 
is suggested, and an anomaly accuracy of 4 to 7 mgal for 1 degree regions has been 
reported (Lerch et al, 1983). 


Table 1 

Goddard Earth Models (Lerch et al 1981 1983) 



Table 2 

Accuracy assessment of ocean geoid (Lerch et al 1981) 


Models 



Comparisons 

GEM-9 

GEM-10 

GEM-1 0B 

GEM-10C 


Geoid heights (m) 






GEOS-3 altimetry - trench areas* 

2.90 

2.87 

2.47 

1.22 

u 

G EOS-3 altimetry - nontrench areas* 

1.92 

1.80 

0.94 

0.75 


Skylab altimetry 

3.2 

3.0 

2.3 



Seasat altimetry 



1.0 



GEOS-3 radial position from 
crossing analysis 

~1 

1.34 

1.00 



Anomalies (mGal) 

G EOS-3 altimeter, 5° blocks 4.7 

Terrestrial, 5° blocks 9.4 

— i I i 

•After bias and tilt have been fit to remove orbit errors. 
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4. NEAR TERM IMPROVEMENTS 


A workshop was held in February 1982 to address future directions for developing 
improved models of the geopotential (NASA, 1982). One recommendation was to develop 
improved models that did not require additional spacecraft to be launched. The con- 
sensus was that perhaps up to a 50 percent reduction in errors up to degree TO, i.e. , 
half-wavelengths of 2000 km, with lesser improvements at higher degree could be 
achieved. These improvements would be realized through two activities. This first 
would be to develop improved modeling tools, i.e.- software incorporating improved 
physical, mathematical and statistical models taking advantage of the enormous compu- 
tational speeds and storage now available. The second is to improve the quality and 
quantity of the existing data by reprocessing, to add additional data from satellites 
not previously used, to incorporate the SEASAT altimetry data, and to collect addi- 
tional laser and radio frequency doppler observations in new campaigns. 

This effort has not yet been formally undertaken. However, even if all expec- 
tations were realized, it would not satisfy the requirements for either orbit deter- 
mination or the ocean geoid to use altimetry for the time-invariant surface 
geostrophic currents. 


5. GEOPOTENTIAL RESEARCH MISSION (GRM) 

Significant improvements in the global geopotential will be possible as a result 
of the GRM, formerly GRAVSAT, which is under study by NASA (Pisacane et al , 1982; 
and Keating, 1983). Accuracy of the global geopotential should be adequate for de- 
termining the sea-surface elevations to 10 cm from altimetry data. This program is 
still in the study pnase with the possibility of a new start in 1988 or later and 
launch in 1992 or later. 

Terrestrial tracking of near-earth spacecraft to refine the geopotential is 
limited. Uncertainties in the propagation velocity, in the ionosphere and especially 
the atmosphere, can induce data reduction errors larger than the orbital perturbations 
of interest. To increase the magnitude of the orbital perturbation It Is necessary to 
decrease the altitude as low as possible. At the 160 km altitude of the GRM space- 
craft, it would require about 276 stations uniformly distributed to provide global 
coverage. At this low altitude, the drag force uncertainty Is about three orders of 
magnitude larger than the gravity forces of interest. These limitations are overcome 
In the GRM by satellite-to-satellite tracking between two satellites in near circular 
polar orbits separated by distances of 100 to 300 km. The disturbing effects of drag 
and radiation pressure and in addition, orbit altitude maintenance can be accomplished 
by the Disturbance Compensation System (DISCOS). This device which was successfully 
demonstrated on an advanced navigation satellite In 1972, TRIAD, uses a mass expul- 
sion system to force the spacecraft to follow the motion of a free proof mass in a 
cavity which Is shielded from the atmosphere and solar radiation. 

Studies have defined the system characteristics given In Table 3. A mission of 
about 7 months, 6 of which will be operational, will require about 1400 ka of hydra- 
zine fuel for each spacecraft, just under half the total mass. An artist's concept- 
ion Is shown In Figure 3. The two spacecraft will be launched by a single shuttle 
mission from the Western Launch Facility. Range-rate measurements between the two 
proof masses will be made to 1 ym^s" 1 at 4 s Intervals. This will be accomplished 
In part by a radio frequency satelllte-to-satelllte doppler system. A laboratory 
Instrument has been developed and tested. Differences with measurements from a 
laser Interferometer were significantly less than the goal, I.e., about 0.03 ynvs" 1 
rms. With a measurement precision of 1 ym.s" 1 at 160 kms altitude, the global geoid 
as a function of degree should be recovered to the accuracy depicted In Figure 4. 



Table 3 
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GRM spacecraft characteristics 


Orbit 

Altitude 

Inclination 

Lifetime 

Launch 

Separation 

Physical characteristics 
Length 

Body diameter 
Solar panels (2) 

Mass 

Power 

Solar panels (2) and body 
mounted cells 

Tracking systems 

Proof mass to proof mass 
Ground based 


~ 160 km 

90 ± 1 deg 

1/2 year (operational) 

Shuttle 

Variable 100 to 300 km 

4.8 m 
0.9 m 

1.5 m x 3.5 m 

2900 kg (1400 kg hydrazine) 

400 W (average) 

1 ium • s' 1 
100 m 





Fig. 3 Artist concept of the Geopotential Research 
Minion. 


Wavelength 
lower bound 
(km) 
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Fig. 4 Geopotential Research Minion (GRM) geoid height 
uncertainty. 



Total geoid height error should be less than 19 cm and satellite altitude error at 
850 km should be about one order of magnitude less. 

The GRM should be able to provide a global geopotential with sufficient accuracy 
to limit the geopotential modeling error of the geoid to less than 10 cm. 


6. SATELLITE G RADIOMETERS 

Satellite-borne gravity gradiometers have been proposed to effect further improve- 
ments in determining the geopotential. This approach is currently envisaged as a 
follow-on to the GRM. To satisfy the objectives of the GRM, measurement precision of 
about 5 x 10" 3 E (IE = 1 0“ 9 s“ 2 ) is necessary at the same spacecraft altitude, 160 km. 

As a result, a measurement accuracy of 10" 4 E over a 4 to 8 s interval is necessary 
to provide an improvement. Reviews of the current state-of-the-art of gravity gradi- 
ometry are given by Wilcox and Scheibe (1983) and Pisacane (1983). 

Current mobile gravity gradiometers can measure to about IE. fundamental limita- 
tions are the instability of the materials, thermal distortion, stability of the 
scale factor and Brownian noise. To achieve the required accuracy it will be neces- 
sary to take advantage of cryogenic technology which should reduce each of the errors 
described above and most significantly the Brownian contribution. An instrument under 
development at the University of Maryland uses two opposed superconducting proof 
masses on a soft suspension and two superconducting sensing coils in a pancake shape 
(Paik, 1981). Two SQUID amplifiers are used to detect the motion from which the 
gravity gradient is determined. A single-axis instrument has been tested and a 
three-axis vector gradiometer should be completed in the near term. A design using a 
superconducting cavity oscillator accelerometer also appears to have promise 
(Reinhardt et al, 1982). Other instruments have been proposed as an interim step to 
achieve 10 -2 E. These are the Bell Aerospace Miniature Electrically Suspended Accel- 
erometer (MESA) and the ONERA proposed CACTUS instrument. 

A spaceborne gravity gradiometer mission in the late 1990's is a possibility. 

The specific improvement in the geopotential will depend on the extent of the 
reduction of errors below 5 x 10“ 3 E. 


7. CONCLUSION 

Spaceborne microwave altimetry has the potential for measuring sea surface 
elevations to a few centimeters. To properly interpret these data for surface 
geostrophic currents, it is necessary to make significant advances in modeling 
the geopotential. The NASA GRM has the potential for meeting the 10 cm accuracy 
requirement for the global geoid. The success of this mission will be a signifi- 
cant achievement in the three-dimensional determination of the geostrophic currents. 
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